Condensed phase of converting boson-fermion mixtures 
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Theory of a condensed state of hybridised bosons and fermions is developed. Normal and anomalous 
Green's functions are obtained diagrammatically and analytically using the Hamiltonian of the 
boson-fermion model (BFM). A pairing of bosons analogous to the Cooper pairing of fermions is 
found. There are three coupled condensates in the model, described by the off-diagonal single- 
particle boson, pair-fermion and pair-boson fields. The Gor'kov expansion in the strength of the 
order parameter near the transition yields no linear homogeneous term in the Ginzburg-Landau 
equation, and the infinite Levanyuk-Ginzburg parameter, Gi = oo, which indicates that previous 
mean- field discussions of BFM are flawed. 
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A two component model of negative U centers cou- 
pled with the Fermi sea of itinerant fermions was origi- 
nally employed to study superconductivity in disordered 
metal-semiconductor alloys 0, 0. Later on it was ap- 
plied more generally to describe pairing electron pro- 
cesses with localisation-delocalisation Q, and to the 
polaron-bipolaron crossover problem in the intermedi- 
ate electron-phonon coupling regime The model at- 
tracted more attention in connection with exotic |5||, and 
high-temperature superconductors |E S Hi IS Oi 
III 111 Q E III 111 111- When the attractive poten- 
tial U is large, the model is reduced to localised bosons 
spontaneously decaying into itinerant electrons and vice 
versa, different from a non-converting mixture of mobile 
charged bosons and fermions [l^ l20l| . More recently this 
boson-fermion model has been adopted for a description 
of superfluidity of atomic fermions scattered into bound 
(molecular) states |2]| . 

Most studies of BFM below its transition into a low- 
temperature condensed phase applied a mean-field ap- 
proximation (MFA), replacing zero-momentum boson op- 
erators by c-numbers and neglecting the boson self- 
energy in the density sum rule |i,|llllil|lll|ll[T8,'2T|. 
When the bare boson energy is well above the chemical 
potential, the BCS ground state was found with bosons 
being only virtually excited 0, ■ MFA led to a conclu- 
sion that BFM exhibits features compatible with BCS 
characteristics ^^^d describes a crossover from the 
BCS-like to local pair behaviour The transition was 
found more mean-field-like than the usual Bose conden- 
sation, i.e. characterized by a relatively small value of Gi 

El 

At the same time our study of BFM |l3 beyond MFA 
revealed a crucial effect of the boson self-energy on the 
normal state boson spectral function and the transition 
temperature Tc- Ref.p^ proved that the Cooper pairing 
of fermions via virtual bosonic states is impossible in any- 
dimensional BFM. ft occurs only simultaneously with the 
Bose-Einstein condensation of real bosons in 3D BFM 
[TsL IT^ . and vanishes in 2D BFM due to the absence of 
the Bose-Einstein condensation in two dimensions |l3l |. 
The origin of this simultaneous condensation lies in a 
softening of the boson mode a.t T = Tc caused by its hy- 



bridization with fermions. The energy of zero-momentum 
bosons is renormalized down to zero at T = Tc, no matter 
how weak the boson-fermion coupling and how large the 
bare boson energy are. Bosons look like overdamped dif- 
fusive modes, rather than quasiparticles in the long-wave 
limit 10, 0, contrary to the conclusion of Ref.Q 
that there is 'the onset of coherent free-particle-like mo- 
tion of the bosons' in this limit. One can expect that the 
boson self-energy should qualitatively modify the whole 
condensed phase of 3D BFM below Tc- 

In this Letter a closed set of equations for fermion and 
boson Green's functions (GFs) is derived taking into ac- 
count self-energy effects in the condensed state of 3D 
BFM. There exist a boson pair condensate along with 
the fermion Cooper pair and the single-particle boson 
condensate in the model. Remarkably, the Gor'kov ex- 
pansion of GFs in the strength of the order parameter 
yields a zero linear term at any temperature below Tc, 
and Gi = 00. It shows that the transition is not a mean- 
field second order transition, and there is no crossover 
from the BCS-like to a local pair behaviour at any values 
of the parameters of BFM. 

The Hamiltonian of BFM in an external magnetic field 
B = V X A is defined as 



H = /'dr^^t(i.)/j(r)V;,(r)+g[0(r)^j(r)^|(r)+if.c.] 

s 

+ i?o0t(r)0(r), (1) 



where V's(r) and 0(r) are fermionic and bosonic fields, 
s =t, J. is the spin, £^0 is the bare energy of bosons with 
respect to their chemical potential 2^, /i(r) = — [V + 
«eA(r)]^/ (2m) — /i is the fermion kinetic energy operator 
, and g is the hybridization interaction converting a boson 
into two fermions and vice versa. Here and further I take 
?i = c = fcs = 1, and the volume of the system V = \. 

The Matsubara field operators, Q = 
exp(i?r)Q(r)exp(-i?T),g = exp(iJT)gt(r) exp(-iJT) 
(Q = ips^'P) evolve with the imaginary time 
-1/T ^ r 1/T as 
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dr 

dr 

d(t>{r,T) 
dr 



hir)ib^{r,T)+g(f>ir,T)4^^{r,T), (2) 

^*(i")V'i(r,T) - .9(^(r,T)V'T(r,T), (3) 
Eo(t)ir,T) + g^l{r,T)^^{r,T). (4) 



The theory of the condensed state can be formu- 
lated with the normal and anomalous fcrmion GFs 
|23,_^(r,r>) = ~{TrMr,r)Mr'fi)), •^+(i-,r',r) = 
{Tripi{r,T)tpi{r' , 0)), respectively, where the operation Tt 
performs the time ordering. Fermionic and bosonic fields 
condense simalteneously Following Bogoliubov 
the bosonic condensate is described by separating a large 
matrix element 0o(r) in 0(r,T) as a number, while the 
remaining part (/)(r, r) describes a supracondensate field, 
(/)(r, t) = 00 (r) + 0(r, t). Then using Eq.(4) one obtains 



g</)o(r) = A(r) = -|^^(r,r,0- 



(5) 



where J^(r,r',r) = {Trtlji{r,T)ipi{r' ,0)). The equations 
for GFs are obtained by using Eqs (2,3,4) and the dia- 
grammatic technique ^4*1 in the framework of the non- 
crossing approximation 25j, as shown in Fig.l and Fig. 2. 

An important novel feature of BFM is a pairing 
of supracondensate bosons, generated by their hy- 
bridization with the fermionic condensate, as follows 
from the last diagram in Fig. 2. Hence, one has to 
introduce an anomalous supracondensate boson GF, 
S+(r,r',T) {Tr(l){r,T)4>{r' ,0)) along with a normal bo- 
son GF, V{r,r',T) = -(T^0(r, r)0(r',O)). The diagrams, 
Fig.l and Fig. 2, arc transformed into analytical equations 
for time Fourier-components of the fermion GFs with the 
Matsubara frequencies a; = TrT{2n+l) (n = 0, ±1, ±2, ...) 
as 

[tij - h{r)]g^{v,v') = S{r - r') - A{v)T+{r, r') 

- 9^TJ2J rfxg_^,(x,r)I?„_„,(r,x)g^(x,r') 

- 9^tY,J dx^+,(r,x)i^^+^,(r,x).F+(x,r'), (6) 

[-iio - h* ir)]T+ (r, r') = A* (r)^^ (r, r') 
- 9^tY, J rfxg^'(r,x)I?^-_^(x,r).F+(x,r') 

+ ff'^E / rfx^--.(r,x)SVc.'(r,x)ec.(x,r')(7) 

{in-Eo)Vn{r,r')=5{r-r') 
9^tY, J dx^..'(r,x)go-c.'(r,x)Po(x,r') 



and, 
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FIG. 1: Diagrams for the normal and anomalous fermion GFs. 
Vertexes (dots) correspond to the hybridization interaction g, 
and zig-zag arrows represent the single-particle Bose conden- 
sate (j)o. 



- g^TY^j dx.F^,(r,x)^o-..'(r,x)B+(x,r'), (8) 
{^zn~Eo)B+{v,v') = 

g^TY, I dxj:-^^,(r,x).FV.'(r,x)Po(x,r') 
- ff'T^ ydxg_^,(x,r)g^,_o(x,r)B+(x,r'). (9) 

for the boson GFs with = 27rTn, and i5(r,r',T) = 
(r,^(r,T)0(r',O)) 

These equations can be formally solved in the homo- 
geneous case without the external field, A = 0. Trans- 
forming into the momentum space yields GFs' time-space 
Fourier components as 



and 



S+(q,c.) = 



iu}* + 



|zcD-aP + |A(k,c.)|2- 

A*(k,w) 
|*cD-ekP + |A(k,c.)|2^ 



in* + Eq 
\in-Eo\^ + \Tiq,n)\^'- 

r*(q,r!) 



\in-Eo\^ + \T{q,n)\ 



2 ' 



(10) 

(11) 

(12) 
(13) 



where Cj = uj + iEj^(k,w), f2 = -I- il]ft(q, fi). The 
fermionic order parameter, renormalised with respect to 
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LD and k, so that it can be absorbed in the renormahzed 
fermion band dispersion. Then the fermion normal and 
anomalous GFs, Eqs.(10,ll) look like the familiar GFs 
of the BCS theory, and one can apply the Gor'kov ex- 
pansion [231 in powers of A(r) to describe the condensed 
phase of BFM in the magnetic field near the transition. 
Using Eq.(7) and Eq.(5) one obtains to the terms linear 
in A 



+ < < 



FIG. 2: Diagrams for the supracondensate boson GFs. The 
Cooper-pairing of fermions leads to the Cooper-pair-like bo- 
son condensate, described by the boson anomalous GF, . 



the mean-field A due to the formation of the boson-pair 
condensate, is given by 

A(k,c.)=A-Kg2r^ J ^T+{k-ci,io')B{ci,uj + cu'), 

(14) 

and the boson-pair order parameter, generated by the 
hybridization with the fermion Cooper pairs, is 

T{ci,n) = g^TY, J ^T{-k,^')T{q^k,n-uj'). (15) 

Ltj' 

Hence, there are three coupled condensates in the model 
described by the off-diagonal fields g(j>o, A, and F, rather 
than two, as in MFA. At low temperatures all of them 
have about the same magnitude, as the fermion and bo- 
son self-energies, 

^f{k,u;) = -g'Tj2 J |^a(q-k,-c.')I3(q,c.-c.'), 

(16) 

(17) 

respectively. 

On the other hand, when the temperature is close to 
Tc (i.e. Tc — T <^ Tc), the boson pair condensate is 
weak compared with the single-particle boson and the 
Cooper pair condensates. Since F, Eq.(15) is of the sec- 
ond order in A, F oc A^, the anomalous boson GF can 
be neglected in this temperature range, where A is small. 
The fermion self-energy, Eq.(16) is a regular function of 



A*(r) = 



dxgi"j (x,r)A*(x)(7(")(x,r). (18) 



The spatial variations of the vector potential are small 
near the transition. If A(r) varies slowly, the normal 

state GF, g!^\r,r') differs from the zero-field normal 



state GF, gi°^(r - r') only by a phase '22^ ^i"^(r,r') = 
exp[— jeA(r) • (r — r')]Q!;^\r — r'). Expanding all quan- 
tities near the point x = r in Eq.(18) up to the second 
order in x — r inclusive, one obtains the linearised equa- 
tion for the fermionic order parameter as 



,(n). 



7[V - 2ieA(r)]^A(r) = aA(r), 



where 



a = 1 



Sb(0,0) 



i-^mii, 

Eo T' 



(19) 



(20) 



and 7 « 7C{S)vj,g'^N{0)/{'i8Tr^T'^Eo). Here vp is the 
Fermi velocity, and N{0) is the (rcnormalized) fermion 
density of states at the Fermi level. 

In the framework of MFA one takes the bare boson 
energy in Eq.(20) as a temperature independent parame- 
ter, Eq = g'^N{0) ln{ii/Tc) 0|, or determines it from the 
conservation of the total number of particles ( the den- 
sity sum-rule) neglecting the boson self-energy I]f,(q, fi) 
H in El IP). Then Eq.(19) looks like the conventional 
mean-field Ginzburg-Landau (GL) equation with a neg- 
ative a oc T — Tc, a relatively small fluctuation region Gi 
[lll |. and a finite Hc2{T) E3- a result one concludes 
that the phase transition is almost the conventional BCS- 
like transition, at least at -Eq ^ 0, ini, and BFM, 
Eq.(l) describes the crossover from the BCS-like to local 
pair behaviour by tuning the parameters |18| . This con- 
clusion is incorrect. The main problem with MFA stems 
from the density sum-rule, which determines the chem- 
ical potential of the system and consequently the bare 
boson energy Eq{T) as a function of temperature. 



)=m- dr|0o(r)p. (21) 



Here t = -1-0, rib = n — Uf is the number of bosons, n 
is the total number of particles, and n/ is the number of 
fermions. The term of the sum in Eq.(21) with Qn = 
is given by the integral 



T 



1 



27r3So + Sb(q,0)' 



(22) 
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where Sf,(q, 0) = T,b{0, 0) + q'^/2M* for a small q is cal- 
culated using Eq.(17) with the normal state fermion GF 
[T^ (here M* is a constant). The integral converges, if 
and only if Eq ^ — l]f,(0,0). This exact result means 
that a ^ at any temperature. In fact, this coeffi- 
cient is strictly zero in the Bose-condensed state, because 
fib — —[Eq + S(,(0, 0)] corresponds to the boson chemical 
potential relative to the lower edge of the boson energy 
spectrum. More generally, /Ltf, = corresponds to the ap- 
pearance of the Goldstone-Bogoliubov mode due to a bro- 
ken symmetry below Tc. As a result, the GL coefficient 
a(T), is zero at any temperature below Tc, a{T) = 0, and 
not only at Tc in the exact theory of BFM. On the other 
hand, MFA violates the density sum-rule, predicting the 
wrong negative a(T) below Tc. 

There are a few important physical consequences. 
Since a{T) = 0, the Levanyuk-Ginzburg parameter 
is infinite. It means that the phase transition is never a 
BCS-like second-order phase transition even at large Eq 
and small g. In fact, the transition is driven by the Bose- 
Einstein condensation of real bosons with q = 0, which 
occur due to the complete softening of their spectrum 
at Tc in 3D BFM. Remarkably, the conventional upper 
critical field, determined as the field, where a non-trivial 



solution of the linearised Gor'kov equation (19) occurs, is 
zero in BFM, Hc2{T) = 0. It is not a finite Hc2{T) found 
in Ref. jl3| using MFA. In the homogeneous case A(T) 
should be determined from Eq.(21) rather than from the 
BCS-like equation (5), which is actually an identity [l^ . 
since Eq = — l]f,(0, 0) below Tc- To get an insight into 
the magnetic properties of the condensed phase one has 
to solve Eqs.(6-9, 14-17) and Eq.(21) keeping the non- 
linear terms. Even at temperatures well below Tc the 
condensed state is fundamentally different from the MFA 
ground state, because of the pairing of bosons. The latter 
is similar to the Cooper-like pairing of supracondensate 
"^He atoms [13, proposed as an explanation of the small 
density of the single-particle Bose condensate in super- 
fluid IIclium-4. The pair-boson condensate should sig- 
nificantly modify the thermodynamic properties of the 
condensed BFM compared with the MFA predictions. 
The common wisdom that at weak coupling the boson- 
fermion model is adequately described by the BCS the- 
ory, is therefore negated by our theory. 
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